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Abstract 

A simple modification of the "Wang-Landau sampling" algorithm re- 
moves the systematic error that occurs at the boundary of the range of 
energy over which the random walk takes place in the original algorithm. 

In two recent papers E] an efficient Monte Carlo procedure was introduced 
that used a random walk in energy space to obtain an accurate estimate of the 
energy density of states g(E) for classical statistical models. If this method 
(now commonly termed "Wang-Landau Sampling") is applied to a restricted 
energy range, effects at the boundaries of the energy range come into play, 
and systematically larger errors in g{E) at the edges of the sampled energy 
interval are observed. Since the method is of quite general applicability, a better 
understanding of these "edge" effects could be of considerable value. Here, we 
show how such an enhancement of errors at the edges can be avoided by a simple 
modification of the algorithm. 

In Wang-Landau Sampling one accepts trial configurations with probability 
min(l, g(E)/g(E')), where g (E) is the energy density of states (DOS) and E and 
E' are the energies of the current and the proposed configuration, respectively. 
At each spin flip trial the DOS is modified g(E) — > g(E) ■ / by means of a 
modification factor /, which is systematically reduced according to / — > f 1 / 2 
whenever the recorded energy histogram H(E) becomes sufficiently flat that 
all entries are within some percentage e of the average energy histogram, i.e., 
H(E) > e(H(E')) e> for all E. H{E) is then reset to zero, and the procedure is 
repeated until a flat H(E) is achieved using a final modification factor f final ■ 
Restricting now the random walk to some sub- interval of the entire energy range 
of the system, one has obviously two basic choices to proceed in case the random 
walk is at the border of the considered energy interval and a spin flip trial would 
result in an energy outside the specified energy segment: 

1. Reject the suggested spin flip and do not update g(E) and the energy 
histogram H{E) of the current energy level E 

2. Reject the suggested spin flip and count the current energy level once more, 
i.e., update g(E) and H(E): g(E) -> g(E) ■ f and H{E) -> H(E) + 1 




Figure 1: Relative error e(E) in g(E) for the first 25 energy levels of a two- 
dimensional nearest neighbor Ising model with linear dimension L = 32. Note 
that the energy scale was not normalized by the number of spins. g(E) was 
obtained by normalizing with respect to the groundstate, e(E) is an average 
over 30 runs. We have used f final — log 10 (8.09 • 10 -10 ) and e = 0.95. 



Method ^ was used in Refs. and this led to a systematic underestimation 

of g(E) at borders of energy intervals 1 . This effect was examined for the two- 
dimensional Ising model, with linear dimension L = 32 and for three different 
ranges of allowed energies: E/(JN) G [-1.7,-1.2], E/(JN) G [-1.8,-1.1], as 
well as E/(JN) G [—1.9, —1.0]. The results showed that systematic deviations 
from the exact DOS occurred only at the right edges of the energy intervals 0. 
Since for the model at hand, this effect only influenced two energy levels directly 
at the border, the recipe used was to overlap the individual intervals over which 
g{E) was sampled by a sufficient number of energy levels so that the affected 
energy levels could be discarded from each when joining the DOS afterwards. 
The asymmetry of this effect can be explained quite simply: For the chosen 
intervals, g(E) has its minimum at the left edge and increases monotonically as 
E approaches the right edge. Hence, during the simulation the random walk is 
"pushed" against the right edge of the sampled energy range, simply because 
generating configurations with energies higher than the right edge energy is more 
likely than generating configurations with energies lower than the boundary 
energy at left edges. Therefore, for each interval, a pronounced effect was only 
visible at the right edge. In order to demonstrate this we have calculated g(E) 
for the first 25 levels of a L = 32 two-dimensional Ising model using single- 
spin-flip Wang-Landau sampling in both variants (method [T] and [2*)l. as well 
as N-fold way updates, which are known not to produce an enhancement of 
errors at edges. In [3] the latter algorithm was tested concerning its behavior at 
edges against the original single-spin-flip version, whereby it was misleadingly 
stated that boundary effects occur when g(E) at edges is sampled the same 
way as inside the energy interval (method [21 . This is actually incorrect. From 
the simulation results, depicted in Fig. one clearly sees that method ^ 
which is almost identical to the implementation of Wang-Landau ^12) leads 

In addition, the energy levels E m i n and E ma x of a specified interval [E m i n , E max ] have 
been updated twice in Refs. when visited by the random walk. 
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Figure 2: Relative error e(i?) in g(E) for the interval E/J e [-80,80] of a 
two-dimensional nearest neighbor Ising model with linear dimension L = 32. 
Note that the energy scale is not normalized by the number of spins. g(E) was 
obtained by normalizing with respect to the exact DOS at the left edge (E/ J = 
—80). s(E) is an average over 5 runs. We have used f final — log 10 (8.09 • 10~ 10 ) 
and e = 0.95. 



to systematic errors in the density of states at the right edge (indicated by a 
dashed line), where two levels are affected, as described in 0. When g{E) is 
sampled according to method [21 no systematic errors are present. In case the 
chosen interval is symmetric around E = 0, the effect should have the same 
magnitude for both edges of the interval, as can be seen from Fig. [5J Again, no 
systematical enhancement of errors at edges is present when the DOS is sampled 
using method 
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